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An analytical description of non-equilibrium phenomena in interacting quantum systems is rarely 
possible. Here we present one example where such a description can be achieved, namely the fer- 
romagnetic Kondo model. In equilibrium, this model is tractable via perturbative renormalization- 
group techniques. We employ a recently developed extension of the flow-equation method to calcu- 
late the non- equilibrium decay of the local magnetization at zero temperature. The flow equations 
admit analytical solutions which become exact at short and long times, in the latter case revealing 
that the system always retains a memory of its initial state. 



I. INTRODUCTION 

Stimulated by the experimental developments in the 
fields of ultracold gases and of nanosystems like quantum 
dots, non-equilibrium phenomena in interacting many- 
particle systems are a fascinating field of current research. 
In the theoretical treatment, a number of conceptual and 
methodological problems arise. Some are related to the 
facts that many interesting phenomena are not perturba- 
tively accessible even in equilibrium, and that moreover 
the limit of a perturbation being small does not neces- 
sarily commute with the long-time limit. Consequently, 
the majority of works on non-equilibrium physics in in- 
teracting systems rely on numerical methods. However, 
those come with their own complications, and analyti- 
cally tractable test cases are desirable. 

It is the purpose of this paper to provide one such test 
case: We shall study the non-equilibrium spin dynam- 
ics in the ferromagnetic Kondo model, described by the 
Hamiltonian 
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in standard notation, with S a and s a (0) being the com- 
ponents of the impurity spin and the conduction electron 
spin density at the impurity site, and j", J 1 - the ex- 
change couplings, with \J ± \ < — J". This model has the 
advantage that its equilibrium physics is well understood 
and can be accessed by perturbative methods, although 
bare perturbation theory needs to be re-summed using 
renormalization group (RG) techniques. Here we shall be 
interested in the decay of a spin prepared in a pure state, 
say | f), after the coupling to the bath is switched on at 
a time t — Q. To make controlled calculations, we shall fo- 
cus on the case of small Kondo couplings \J"\, \ J ± \ <D 
where 2D is the conduction-electron bandwidth. While 
the equilibrium physics of Eq. ([T]) is that of an asymptot- 
ically free spin, the non-equilibrium problem is far from 
trivial. In particular, it is not obvious whether the expec- 
tation value {S z )(t) relaxes to zero: The equilibrium flow 



of the spin-flip coupling towards zero is slow (i.e. loga- 
rithmic), which suggests that spin- flip events occur even 
at long times. In case the spin relaxes to a non-zero value, 
the follow-up question is whether this value is identical 
to the equilibrium expectation value (S z ) cq obtained in 
the presence of an infinitesimal magnetic field. 

In order to tackle this problem, we shall employ a non- 
equilibrium extension of the flow-equation method, orig- 
inally developed by Wegner— and independently - in the 
field of high energy physics - by Glazek and Wilson. 13 
For equilibrium settings, this method utilizes continu- 
ous unitary transformations to successively eliminate off- 
diagonal terms in the Hamiltonian, starting from the 
highest energies. Although somewhat similar in spirit 
to RG techniques, the flow-equation method has the 
advantage of diagonalizing, instead of eliminating, the 
high-energy part of the Hamiltonian, such that the en- 
tire system is described at any stage of the transforma- 
tion. The flow-equation method can be efficiently em- 
ployed to describe the non-equilibrium time evolution of 
observables^ The Hamiltonian is unitarily transformed 
into a diagonal form, together with the observables of 
interest. In this representation, the time evolution can 
be determined exactly. Finally, the result is transformed 
back using the inverse of the flow-equation transforma- 
tion. In this way, the accumulation of errors in the long- 
time limit is prohibited, and secular terms do not occur. 
In Refs. [HIH, this method was applied to dissipative quan- 
tum impurity models, namely the dissipative harmonic 
oscillator and the ohmic spin-boson model, and excellent 
agreement with available analytical and numerical results 
was obtained. 



A. Summary of results 

In the following, we summarize our main results ob- 
tained in the limit of weak coupling, distinguishing the 
cases of isotropic and anisotropic Kondo couplings. Most 
importantly, in both cases the relaxation of (S z (t)) 
as t — > oo is neither to zero nor to the equilibrium 
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value {S z ) eq . Instead, the asymptotic value is given by 
1/2- (S*(t -» oo)) = 2(1/2- (S z } oq ), i.e., the spin polar- 
ization is reduced twice as much compared to the equi- 
librium case. 

The initial (short-time) decay of (S z (t)) is set by J- 1 
(independent of jH), (S z (t)} = ± [l - (2J ± Dt) 2 ] . 6 In the 
isotropic case, the long-time decay towards the asymp- 
totic value is logarithmic, (S z (t)) — (S z (t — > oo)} oc 
1/lnt. In contrast, for anisotropic couplings with J" < 
— | J | the logarithms in the long-time limit are replaced 
by power laws, {S z (t)) — (S z (t — > oo)} oc i 2 -'", where 
j II < is the dimensionless fixed-point value of j" . These 
results do not come unexpected, considerung the equi- 
librium RG flow for the ferromagnetic Kondo problem, 
which is characterized by logarithmic (power-law) flow in 
the isotropic (anisotropic) case. 

Below, we shall give explicit expressions for (S z (t)) for 
all times, Eqs. (|42|) and (|48p. and discuss in detail their 
evaluation in both the short-time and long-time limits. 



B. Relation to earlier work 

Our work adresses both recent methodological devel- 
opments and fundamental theoretical questions. Ap- 
plications of the flow-equation approach to interaction 
quenches in quantum many-body systems have been ap- 
plied in several recent works studying dissipative quan- 
tum systems 2 -^ and interaction quenches in the Hubbard 
model. 7 '- 8 These works relied on either numerical imple- 
mentations or analytical approximations restricted to fi- 
nite time scales. Previous numerical results exist for 
the real-time dynamics of the antiferromagnetic Kondo 
model, which has also been treated with the flow equa- 
tion method^. The ferromagnetic Kondo model has been 
discussed numerically with the time-dependent Numeri- 
cal Renormalization Group (TD-NRG) method^ at finite 
temperatures^ or in the context of the underscreened 
Kondo effect^ These works did not identify the asymp- 
totic long-time behaviour of this problem at T — 0. This 
behaviour has been adressed by a comparison of the flow 
equation approach and the TD-NRG in a recent letter 4 
The asymptotic long-time tails have been clarified by the 
analytical asymptotics of the flow-equation approach, to 
be described in more detail in the present paper, and the 
numerical implementation of the TD-NRG. 

Relaxation of an impurity spin in a host model was 
already discussed by Langreth and Wiikinsjii who de- 
veloped a theory of spin resonance in dilute magnetic 
alloys. Based on the Kadanoff-Baym approach, these 
authors derived Bloch-like equations for the paramag- 
netic resonance, which predict a relaxation of the dis- 
turbed magnetizations to the instantaneous local equi- 
librium magnetization. Our work differs in predicting 
a relaxation towards a local non-equilibrium magnetiza- 
tion under slightly different conditions than considered 
by Langreth and Wilkins. A perturbative expansion 
of Kadanoff-Baym equations as used by Langreth and 



Wilkins would fail to reproduce our result, since bare per- 
turbation theory cannot reproduce our asymptotic long- 
time relaxation laws. 

C. Plan of the paper 

The remainder of the paper is organized as follows: 
In Sec. |TT] we start our considerations by studying an 
exactly solvable toy model of a spin coupled to two 
fermions. This model illustrates a number of remark- 
able features which will re-appear in the treatment of 
the Kondo model. Sec. IIIII describes the flow-equation 
transformation for the ferromagnetic Kondo model, in- 
cluding the transformation of the relevant observables. 
The explicit calculation of the time-dependent impurity 
magnetization is subject of Sec. IIVI The central results, 
obtained both analytically for short and long times and 
numerically from a full solution of the flow equations, are 
presented in Sec.fV] A brief discussion of limitations and 
applications closes the paper. Technical details will be 
relegated to the appendices. 

A brief account of our results has been published in 
Ref. H together with a comparison to numerical data 
obtained by the TD-NRG method. 

II. TOY MODEL 

Several interesting features of the non-equilibrium pro- 
cesses governed by the ferromagnetic Kondo model can 
be captured by a simple exactly solvable toy model con- 
sisting of two fermion levels coupled to a spin 1/2 impu- 
rity spin 5* through an SU (2)-symmetric exchange cou- 
pling: 

H = J2 ( C « C « - d M + f I>« + d aK/3( c /3 + d p) ■ 
a a,/3 

(2) 

Due to the finite Hilbert space the real time evolution 
problem for an arbitrary initial state becomes exactly 
solvable. Specifically we look at \ip(t = 0)) = \ipo), where 

IVo) d = 10} ® | u> ® I T) (3) 

denotes an unoccupied c-electron level state |0), a dou- 
bly occupied d-electron level state | TO and the spin-up 
state | t) of the impurity spin. Notice that \4>o) is the 
ground state of { c a c a ~ d a^a) ~ B S z , where B > 
shall cause an infinitesimal Zeeman splitting. Later we 
will generalize the initial state to be a product state of a 
Fermi sea with the impurity spin-up state for the dynam- 
ics generated by the actual ferromagnetic Kondo model. 

In our toy model we now investigate how the spin ex- 
pectation value deviates from its initial value at time 
t = due to the dynamics generated by @, that is we 
study the observable 

O = S z \ . (4) 
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Let us denote the exact eigenstates of H with eigenen- 
ergies E n by \ip n )- Their explicit construction is shown 
in appendix El in fact one just needs to diagonalize a 
3 x 3-matrix. 

Since the initial state l^o) is not an exact eigenstate, 
the expectation value of O acquires a time dependence 
given by 



0(t) = ^ \e Mt Oe- lHt \^ ) 



U n Un ' 



i(E n ,-E n )t 



where we have used the decomposition \ipo) = 
^2 n u n \ip n ) . For an examplary case 0(t) is depicted in 
Fig. [1] one can clearly observe oscillations with a finite 
number of Bohr frequencies E n > — E n . Since there are no 
accidental degeneracies, the time average is given by 



dt(il> \e im Oe- im Wo) 



0(t) = lim 



/ x T 

|2 / 



(5) 



We compare this average with the ground state expecta- 
tion value in the interacting system ([2]): 



Oeq = (ljjeq\0\lp eq ) , 



(6) 



where \ipeq) is the exact eigenstate with the smallest 
eigenenergy E n . The ratio 0(t)/O eq is depicted in Fig. [2] 
as a function of the coupling strength. We observe that 
generically the ratio differs from one, which shows that 
the time-evolved initial state \ipo(t)) always retains a 
memory of its initial preparation. While this is not at all 
surprising for a finite Hilbert space like in our toy model, 
we will later see that this observation is generalized to 
the ferromagnetic Kondo model in the thermodynamic 
limit. 

The second observation in Fig.[2]is the universal value 2 
in the weak-coupling limit. The universality of this result 
is guaranteed by a theorem proven in Ref. 8: Take a dis- 
crete quantum system governed by H = Hq + g Hi nt , 
where g Hi nt is a weak perturbation for which non- 
degenerate perturbation theory is possible. Let |0) be the 
ground state of H and |0) be the ground state of H. Let 
O be an observable that commutes with Hq and which 
annihilates |0). Then 



a. (Ol^Oe^lQ) 
(0|O|0) 

One can easily verify that the conditions of this theorem 
are met in our toy model if we take 

Ho = 52{cU a -dUa)-BS* 

a 

Hint = ^-^(4+4)^/3(0/3 + ^) 

a, j3 
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FIG. 1: The time dependence of the expectation value 
(if>o\S z (t) - | |^o> is shown for g = -1/2. 




<V. |S*(t)-l/2|V> > 



FIG. 2: We depict the ratio r(g) — {s^-i/2) e 

function of the coupling strength g. In the weak-coupling 
limit one finds the universal ratio r(0) = 2, see text. 



with infinitesimal B > 0. In appendix [X] we show ex- 
plicitly how this theorem comes about in our simple toy 
model. 

Again we will find that this universal factor 2 between 
the time-averaged expectation value in the time-evolved 
initial state and the equilibrium ground state also holds 
in the weak-coupling limit of the ferromagnetic Kondo 
model. 



III. KONDO MODEL AND FLOW-EQUATION 
TRANSFORMATION 



In this section, we summarize the equilibrium proper- 
ties of the ferromagnetic Kondo model which are relevant 
for the subsequent discussion. We then explain in some 
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detail the flow-equation treatment of the model, which 
shall be used in Sees. IIVI and [V] to calculate the non- 
equilibrium magnetization. 



A. Definition of the model 

The Kondo Hamiltonian, Eq. |T|), can be re- written as 
H = Yercl c t 

*- — / fc ka ka 

ka 

Z-~t k'k k'k 

(8) 



k,k' 



k,k' 

Here, S ± = S x ± iS v , and the conduction electron spin 
density is given by 



S k'k 



9 k'a 



a, /3 



z,± 



(9) 



where <r ± = \{ci x ± ia y ) derive from the Pauli matri- 
ces. The anisotropic exchange couplings Jjk~ and jjL 
usually describe scattering on an isotropic Fermi surface, 
and their momentum dependence can then be safely ne- 
glected. In physical results, they enter proportional to 
the density of states at the Fermi surface (pf), thus 

we define the dimensionless couplings j" == p^j" and 
3 = PfJ^- 



B. Low-energy fixed points 

The important difference between the ferromagnetic 
and the antiferromagnetic regime of the Kondo model 
can be understood from Anderson's poor man's scaling 
analysis. 1 In this approach, the band width of the con- 
duction electrons is progressively reduced by a flowing 
cutoff energy A and the renormalizcd interactions due to 
the elimination of virtual excitations to the band edges 
are calculated perturbatively. As shown by Anderson, 
this procedure leads to the scaling equations 



dhxA 
dJ 1 - 
rflnA 



- PF J ±2 
-PfJ x J ]1 , 



(10) 



where the scaling trajectory for ferromagnetic couplings 
with \J ± \ < — J" is given by the conserved quantity j" 2 — 
J ±2 — const. 

For antiferromagnetic couplings, more precisely for 
jll > or |J X | > — jH, the RG flow is driven to strong 
coupling, i.e., j" — > oo and jj^j — > oo, and the scaling 



equations have no low-energy fixed point at finite cou- 
pling. On the level of perturbation theory, logarithmic 
divergencies ln(fcsT/D) occur at low temperatures. 

In contrast, in the ferromagnetic Kondo model with 
| < — J", the perturbative renormalization of the cou- 
pling constant remains controlled in the limit A <C T 
since the coupling monotonously renormalizes to a fi- 
nite value J(T) in the limit A — > 0. At T = and for 
anisotropic couplings, the longitudinal coupling remains 
finite and the transverse coupling renormalizes to zero 

according to the power-law J X (A) oc A PF V^ . In the 
isotropic case, both couplings logarithmically renormal- 
ize to zero following J (A) = J/(l + pfJ^(A/D)), and 
the impurity becomes asymptotically free at low ener- 
gies. This important property makes it possible to use 
the flow equation renormalization scheme in a perturba- 
tively controlled manner in the following. 



C. Flow-equation method 

Let us briefly review the basic ideas of the flow- 
equation approach (for more details see Ref. [H). A 
many-body Hamiltonian H is diagonalizcd through a se- 
quence of infinitesimal unitary transformations with an 
anti-hermitean generator fy(-B), 



dH{B) 
dB 



= [ V (B),H(B)} , 



(11) 



with H(B = 0) the initial Hamiltonian. The "canon- 
ical" generator— is the commutator of the diagonal 
part Hq with the interaction part iJj nt of the Hamilto- 
nian, rj(B) = [Hq(B), fli n t (S)]. Under rather general 
conditions the choice of the canonical generator leads 
to an increasingly energy-diagonal Hamiltonian H(B), 
where interaction matrix elements with energy trans- 
fer AE decay like exp(— B AE 2 ). For B — > oo the Hamil- 
tonian will be energy-diagonal and we denote parameters 
and operators in this basis by e.g. H = H(B — oo). 

The key problem of the flow equation approach is 
generically the generation of higher order interaction 
terms in Eq. pip , which makes it necessary to truncate 
the scheme in some order of a suitable systematic expan- 
sion parameter (usually the running coupling constant). 
Still, the differential nature of the approach makes it pos- 
sible to deal with a continuum of energy scales and to de- 
scribe non-perturbative effects. This has led to numerous 
applications of the flow-equation method where one uti- 
lizes the fact that the Hilbert space is not truncated as 
opposed to conventional scaling methods. 

In Refs. [US these features have been utilized to 
treat the real-time evolution of observables under non- 
equilibrium conditions. The general setup is described 
by the diagram in Fig.[3l where \ipi) is some non-thermal 
initial state whose time evolution one is interested in. 
However, instead of following its full time evolution it 
is more convenient to study the real time evolution of a 
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B = 



B = oo 



H,A, 



U 



Non-perturbative solution 
of Heisenberg equations 
of motion for operator A(t) 



u + 









Time evolution 


A{t) = e^'Ae-™, |*,) 



FIG. 3: The forward-backward transformation scheme in- 
duces a non-perturbative solution of the Heisenberg equations 
of motion for an operator. U denotes the full unitary trans- 
formation that relates the B = to the B = oo basis.— 



given observable A that one is interested in. This is done 
by transforming the observable into the diagonal basis in 
Fig. [3] (forward transformation) : 



dO 
dB 



[r)(B),0(B)} 



(12) 



with the initial condition 0(B = 0) — A. The central 
observation is that one can now solve the real time evo- 
lution with respect to the energy-diagonal fl exactly, 
thereby avoiding any errors that grow proportional to 
time (i.e., secular terms): this yields A(t). Now since 
the initial quantum state is given in the B = basis, one 
undoes the basis change by integrating Eq. (JTSJ) from 
B — > oo to B — (backward transformation) with the 
initial condition 0(B — > oo) = A(t). One therefore effec- 
tively generates a new non-perturbative scheme for solv- 
ing the Heisenberg equations of motion for an operator, 
A(t) = e iHt A{0)e' lHt . 



D. Flow equations for the Hamiltonian 

It is a straightforward calculation to apply the flow 
equation technique in the outlined manner to the ferro- 
magnetic Kondo Hamiltonian i 12 ! 16 For the sake of sim- 
plicity, we consider the case of S = 1/2 here and give 
a generalization to arbitrary spin S in appendix [Bj To 
start with, we split up the flowing Kondo Hamiltonian as 

fl(fl) =fl + flint (fl) (13) 

where H — YV e,-*ct C - , and H int (B) describes the 

* — * Ho h klJ k(J 

interaction with the spin-1/2 degree of freedom S, 

hub) = E4^ :s k : 

k'k 
k'k 

Normal ordering of the fermion operators has been de- 
noted by : :. It is used in order to safely truncate 
normal-ordered terms with respect to a given refer- 
ence state from the ansatz Eq. (fT^tj) . These truncated 



terms are in addition proportional to higher powers of 
the coupling J compared to the leading order of the 
calculation^ In the following, we will use normal or- 
dering w.r.t. the equilibrium density matrix of the non- 
interacting Fermi sea, but other choices could be used as 
we\\&H According to the ansatz Eq. (fTl| . we will cal- 
culate the flow of the Hamiltonian only to first order in 
the coupling J. It would be straightforward to formulate 
an ansatz that also contains coupling terms proportional 
to higher orders in J, which are necessarily generated 
during the flow. For our purpose, it will be sufficient 
to work only to lowest nontrivial order, which includes 
only interaction terms of first order in the coupling J. 
The canonical generator is immediately obtained from 
[An, #i n t (£?)], 

T)(B) = [H Q , fli n t (A)] 



with 



k'k 

4'k^ 



' S k'k ' 



k'k 



(15) 



nL(fl) d = tea -e : )jljB) 

l k , k \ J \ k' k> k i k \ I 



clef 



By comparing coefficients on both hand sides of the dif- 
ferential equation 4§ = \q{B), fl(fl)], one immediately 
finds the flow equations for the coupling constants^ 



dB 1 k ' k) k'k 

+ E( 2 ^-^-^) J fe Vi 



x (--n(q)) + 0(J 3 ) 



d 4k 



dB - y £ k' £ k ) J k 'k 
+ Y.^ e ^ e k'- £ k) 



x i (A J± + Ji A) (1 - n(q)) 

2 V k'q qk k'q q*k> \ 2 

+ 0(J 3 ) . 



(17) 



The initial conditions for the exchange interactions are 
J l L(B = 0) = JH /N and J±AB = 0) = J ± /N, where N 
is the number of band states. In the flow equations (1 1 7[) . 
the Fermi distribution function n(k) = l/(l+exp(£g/T)) 
has been introduced, which enters the flow equations as a 
contraction generated by our normal-ordering procedure. 
Since the couplings Jgt ^ ow t° zer0 m ^ ne nm it B — ► oo, 
the fixed point of the flow-equation renormalization of fl 
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is 

H(B = 00) = ^ e ^ c L c k<T 

kcr 

+ E J k k (S = ^V c Lc kCT . (18) 

kcr 

For a discussion of the magnetization curve (S z (t)), 
we can neglect the term J2k<r Jkk(B = °°)^ Zac ka c k<7^ 
since the renormalized couplings </J[ k (-B = 00) shift the 
fermionic energy levels proportional to the inverse num- 
ber of band states (1/N) and will drop out of the impurity 
magnetization curve (S z (t)) in the thermodynamic limit 
N — > 00. It is therefore possible to exploit the nonin- 
teracting form H(B = 00) = J^ka e kcrC kcr c k(T by trans- 
forming the impurity spin operator S into the same basis 
representation. 



with a = (x, y, z) and the initial condition S a (B = 0) = 
S a . By considering the commutator [rj(B), S a ], it is 
straightforward to work out the following ansatz for the 
flowing operator S a (B) 

S a (B) = h{B)S a + J2y~ k (B) : (S x sg^ :, (20) 

k',k 

which is accurate to neglected normal-ordered contribu- 
tions of 0(J 2 ). In the following, we transform only the 
operator S z which describes the impurity magnetization. 
By extending the coefficients in the ansatz of Eq. (|20|) 
by a dependence on time, the real time evolution of the 
operator S Z (B) can be parameterized. We will make use 
of this parameterization later on, and therefore introduce 
here the following abbreviations for the couplings, which 
are related to specific values of the parameters B and i: 



E. Flow equations for the spin operator 

In order to transform the impurity spin operator S 
into the basis of the diagonal Hamiltonian H(B = 00), 
we need to solve the differential equation 



dS a {B) 
dB 



[r)(B),S a (B)], 



(19) 



h(B,t) = < 



h(t) B 

h(B) t = 

h(t) B 

h B 



co, t = 



The differential flow of the ansatz Eq. (|20|) is readily 
obtained from the commutator [r)(B), S Z (B)] as^ 



% - ^(e g ,-e u )J^(B) 7w (B)n(k')(l-n(k)) 

k'k 

£jg£ = h(B)(sg edJfoB) - \ 2 {{eg ~ ^J^h^) + (e S - s^B^B)) (1 - 2n(u)) . (21) 



In Sec.fV] these differential equations will be simplified in 
order to derive an approximate solution to them. For a 
check of approximations, it will be useful that the flowing 
couplings h(B) and 7gg, (B) are related by the sum rulei^ 



(S Z (B) 2 ) = 5 (/ l 2 (B) _^ 7 „ mn( fe 0(1 _ n( fc )) ) 

kk' 

= \ + 0{J 2 ) , (22) 



where (S Z (B) 2 ) is evaluated with respect to the non- 
interacting Fermi sea at equilibrium and a polarized im- 
purity spin. This sum rule can be easily proven by dif- 
ferentiating it with respect to B and inserting the flow 
equations (pH) . 



IV. TIME-DEPENDENT MAGNETIZATION 

Since a ferromagnetically coupled Kondo spin is 
asymptotically free at low energies, its equilibrium prop- 
erties are only perturbatively renormalized in comparison 
to a free spin. In the case of isotropic couplings, it has 
been shown by Abrikosov and Migdal 1 ^ by diagrammatic 
means that for general spin quantum number S the equi- 
librium magnetization of the impurity (M = gi/iB(S z )) 
in a magnetic field £?, caused by the Zeeman term 

Hzeeman = _ SiMS'S' B 

is renormalized to one-loop order as 

Here, we introduced x = ln(e p / '{gi^sB)) with the Fermi 
energy ef, the gyromagnetic ratio gi of the impurity and 
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the Bohr magneton fig. Spin- flip scattering effectively 
reduces the magnetization but is not able to screen the 
spin completely. If otherwise an impurity spin is decou- 
pled from the conduction electrons and thus its magne- 
tization is completely polarized, a subsequent coupling 
to the conduction band will lead to spin-flip scattering 
and the impurity magnetization will start to decrease as 
a function of time. In order to analyze this behaviour 
quantitatively, we assume in the following that the com- 
plex of impurity and conduction electrons is prepared in 
the product inital state 



= | T) ® \FS) , 



(24) 



where \FS) is the non- interacting Fermi sea at equilib- 
rium. We assume furthermore that impurity and con- 
duction electrons become coupled instantaneously at the 
time t = 0. 

Using a perturbative expansion of the Heisenberg equa- 
tion of motion for the magnetization operator S 12 , it is 
straightforward to calculate the perturbative time evolu- 
tion of the magnetization as (with (S z ) = \)r^ 

/OO />OG 
de / de'p(e) P (e') 

-I - cos((e - e')t)- 



x f{e'){l-f{e)) 
+ 0(J 3 ), 



e') 2 



(25) 



with the Fermi function /(e) and the conduction elec- 
tron density of states p(e) — J2k <K e ~ e k) /N, where N is 
the number of band states. This result neglects contri- 
butions of 0( J 3 ) in coupling strength. It is a well known 
property of the Kondo model that renormalization of the 
couplings J 1 - by j" sets in precisely at third order in J.— 
At sufficiently large time scales, the low-energy couplings 
are expected to dominate the magnetization dynamics, 
and the renormalization of spin-flip scattering will be- 
come important. Beyond short times, the perturbative 
result of Eq. (|25|) is therefore expected to break down. 

In the limit T — > and for a flat band with density of 
states pp = 1/{2D)^ the integrals can be rewritten as 

(S z (t)) = (S z ) - (S z ) \ W - cos(xt))dx, (26) 

Jo % 

where we introduced the function 



x < D 

K W - <j 2D-x, x > D. 



(27) 



This integral can be solved analytically with the result 

(S z (t)) = (S z ) ~ (S z )(2J ± ) 2 [G(2Dt) - 2G{Dt)] (28) 
where the function G{x) is defined by the series expansion 



G(a;) = |-J (2iy.2i(2i-i) x21 - 



(29) 



As argued above, Eq. (j28|) fails to predict the correct 
saturation behaviour of the magnetization, see also Figs[H 
and [6] 

To improve upon this unsatisfactory result, we perform 
the perturbative solution of the Heisenberg equation for 
the operator S z in the basis given by the diagonalized 
form of the Hamiltonian, see Fig. [3J In several previous 
applications^^ this approach has been used to obtain 
the real-time evolution of physical observables well be- 
yond the perturbative short-time regime. 

In the corresponding new basis representation, the so- 
lution to the Heisenberg equation of motion for S z is 
readily obtained from S z (t) = e lHt S z e~ lHt . Since the 
operator S z is a conserved quantity in this basis, in the 
ansatz Eq. ([20]) only the coefficients 7££, will obtain time- 
dependence, with the result 



Sz(t) 



7pg(<) 



hS z 

E?«(*) 

k',k 

w i4(£ *'- £s) 



(30) 



The inversion of the unitary transformation (denoted by 
U, see Fig. [3]) will yield the effective non-perturbative 
solution of the Heisenberg equation of motion for the 
operator S z , which we write as 

S z (t) = fc(*)S" + EW*) = (^ s + ^i) 

k',k 

+ 0(J 2 ). (31) 



All coefficients in Eq. (|3Tj) are obtained by integrating 
the flow equations (f2Tj) from B — » oo to B = 0, with the 
initial condition in B — > oo posed by the operator S z (t), 
Eq. (001) ■ From Eq. ([3TJ). the magnetization (S z (t)) 
can be readily read off as (S z (t)) = h(t)/2. It remains 
to determine the coefficient function h(t), for which an 
analytical result can be derived in a limit of long times. 



V. ANALYTICAL RESULTS FOR THE 
MAGNETIZATION 

It is too ambitious to solve the flow equations (fT7|) 
and (|2ip analytically as they stand, but much progress 
can be made by choosing a suitable ansatz for the flow of 

the coupling constants jj and ji-, . Several suggestions 
along this line have been given in Ref. Il2|, chapter 2. In 
the limit of exchange couplings J—" close to the Fermi 

kk' 

surface, meaning \e^\/D, \eg\/D <€. 1, the flow of these 
couplings can be formally parameterized as 



(32) 



We call this parameterization of the flow the infrared 
parameterization with the infrared couplings J^\B), 
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since it is asymptotically exact for £z,£a = 0. The 

flow of the couplings J^{B) is identically to that of 

Jfep'l (B), which are the couplings at the Fermi surface.— 
Since during the flow, high energy couplings with en- 
ergy transfer AE are eliminated according to the relation 
AE oc S -1 / 2 , the flow of the low-energy couplings can 
be shown to be negligible for B < D~ 2 J2. By suitably 
adjusting the initial conditions to J"(B = D~ 2 ) = J" 
and J ± (B = D~ 2 ) = J 1 - the flow of the dimensionless 

coupling constants j" = PfJir an( l = PfJjr coin- 
cides then exactly with the scaling equations (fTtJf , where 
the flowing bandwidth A is identified as A = B^ 1 / 2 . 

In addition, it is possible to simplify the flow equa- 
tions (|2ip for the spin operator S z by neglecting the con- 



tribution to the derivatives Z k J> k of second order in the 



juplings J" and J -1 , setting 



dB 



dB 



= h(B)(eg - e % )jUB) + 0((eg - e- k )J 2 ) . (33) 



This approximation has to be carefully justified for 
anisotropic couplings \ J ± \ < |J"| and jH < 0, since the 

couplings J- - contained in the neglected terms flow to a 

k'k 

finite value at the Fermi surface (given by eg = eg=0), 
while Jg-.(B) will flow to zero. We justify this approxi- 
mation in detail in appendix [Cl 

Our approximations so far aim at the behavior of the 
couplings 7£g, near the Fermi surface, which therefore 
are suitable to obtain the long-time asymptotics of the 
impurity magnetization. We will show now that these 
low-energy couplings depend decisively on the behavior 
of the scaling equations (fTO)) , leading to qualitatively dif- 
ferent results in the anisotropic regime vs. the isotropic 
regime of the exchange couplings. 



A. Isotropic regime 

In the following, we restrict ourselves to the isotropic 
case J = J 1 - = jH. 



1. Flow of the dimensionless coupling j 

As already pointed out before, we will be mainly inter- 
ested in the behaviour of the flowing coupling JgAB) for 
energies close to the Fermi surface. Exactly at the Fermi 
surface, the flow of the dimensionless coupling j = Jpp 
according to the scaling equation (fTU|) can be integrated 
exactly, yielding 



m = 



(34) 



We conclude that the flow of the infrared coupling is given 
by Jir(B) = J/(l-(j/2) ln(BD 2 )), using the correspon- 
dence A = B~ 2 . 



2. Couplings of the spin operator 

For a further simplification of the flow equation (|33l) . 
it is important that the coupling h will be a number close 
to 1, with corrections that vanish continuously as J — > 0. 
Using this fact, we set h(B) = 1 in Eq. (|33j) . leading to 
corrections that become small for small J, as we explain 
now. 

More precisely, this approximation can be justified by 
a wellknown approximation for the impurity magnetiza- 
tion in an external magnetic field. In presence of an in- 
finitesimal Zeeman splitting BS Z , B — > + , the coupling 
h describes the ground state expectation value of the op- 
erator S z by the relation h/2 = (S z ). By multiplying 
(S z ) = h/2 with the impurity gyromagnetic ratio gi and 
the Bohr magneton /is, the impurity magnetization is 
given as M = gifiBh/2. However, this magnetization can 
be easily read off from (|23|) by taking the limit B — > + , 
which is M = (l/2)gfj//s5(l + j/2), and thus we obtain 
(for S =1/2) 



(35) 



where possibly corrections of 0{j 2 ) occur. Since h is 
renormalized only to O(J), equation Q33p can be simpli- 
fied by setting h(B) = 1 and neglecting corrections of 
0{J 2 ). 

We now aim at the flow of the couplings ^gt{E), 
which can be evaluated asymptotically in the limits 
|eg|, \eg \ < D and (D - e % )/D « e % ,)/D < 1. 

First, we concentrate on couplings near the Fermi surface 
(led, \eg | <C -D), for which the infrared parameterization 
from Eq. (|32|) can be used in Eq. (|33[) . and we obtain 



= fin 

Jo 

+ 0(J 2 ) 



e g) J -^pe-^-^ d B> 



(36) 

According to Eq. (1341) . the infrared couplings Jm(B) 
depend only logarithmically on B, and in the integral 
©, we can set Jir(B) = J IR (B = (eg - eg) ) + 
0(J 2 ), leading to the result 

kk N{eg-ez) ' 

+ 0{J 2 ) 

Jir{B= (eg-e^y- 



0{J 2 ) 



(37) 



Now, it is possible to determine the coupling h as a so- 
lution of Eq. QTJ. Employing Eq. ([33]) in Eq. ^> and 
setting h(B) = 1, we first obtain the formal result 

1 



h-1 = - 



dB—^giBMk'Kl-nik)) 



2^7 ( 

z - - Jo 

kk' 

\jllh n (k'){l-7l{k)). 

kk' 



(38) 
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At low temperatures, only the low-energy couplings from 
Eq. (|36p contribute, and at zero temperature and using a 
flat band with density of states pf = jjy, we can rewrite 
the momentum sums as an integral. Employing finally 
the logarithmic scaling behavior of the infrared couplings 
J IR (B = (eg - e^y 2 ) from Eq. JTDJ) in Eq. (SSJ), we 
arrive at the result 



2D 



dx 



k{x) 



3 



(l+jln(§)) 



(39) 



This integral can be evaluated exactly, and to leading 
order in j we obtain the result h = 1 + j/2 + 0(j 2 ). 
To leading order in j, this result is confirmed by Eq. 
(|35|) . which has been originally derived by Abrikosov and 
Migdal^ 



3. Time-dependent magnetization 



As discussed in Sec. [TV] the time-dependent magneti- 
zation (S z (t)) follows directly from the coefficient hit). 
It is therefore necessary to solve the flow equations (|2Tj) 
for time-dependent initial conditions, which endows the 
flowing couplings with an additional time-dependence, 
and we change our notation according to h(B) — > h(B, t) 

and7 S p(B)-7 S p(B,*). 

Nevertheless, it is trivial to obtain jrg {B, t) to leading 
order in </, since the derivatives in Eq. (137|) do not depend 
on time and can be integrated in analogy to Eq. ([57]) . 
with the result 

7*s(B,t) = 7Pjf(*)+ / ~ £%)J9%{B>)dB' 

J CO 

+ 0{J 2 ) 

+ 0(J 2 ) . (40) 

This expression can be employed in the flow equation (|2ip 
for the coupling h(B, t), which can be integrated in anal- 
ogy to Eq. (|38|) . The formal result for the coupling h(t) 
is finally 



h(t) 



kk' 



- 1) 



x n(k')(l - n{k)) 



dB 

kk' 



x n(fc')(l -n(fc)) 



(41) 



In appendix [D] we derive this result in a different way by 
employing the sum rule from Eq. (|22|) , thereby further 
justifying all approximations. Aiming at the long-time 
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FIG. 4: We show (S z (t)) for isotropic coupling j and short 
times. For very short times t <C -D -1 , the perturbative re- 
sult from Eq. (|28[) (dashed line) is asymptotically coinciding 
with a fully numerical solution of the fkow equation (I21|l (full 
lines). The renormalization of J x by j" sets in at energy 
scales E < D, leading to a reduction of spin-flip scattering. 
Beyond the short-time regime t <§C D _1 , the magnetization 
relaxes therefore slower than predicted by the perturbative 
result. 



TABLE I: For large times, we fitted the full numerical solu- 
tions of the flow equations against the function 0.5(l + aj + 
l/(lnt+l/(&J))) . For couplings — j < 0.1, quantitative agree- 
ment with the analytical result (S z {t)) = 0.5(l/(lnf- 1/j) + 
1/2+ j/2 + 0{j 2 )) is very , 







-0.05 


-0.1 


-0.15 


a 




0.441 


0.442 


0.441 


b 




0.418 


0.400 


0.382 


behaviour at 


T = 


= 0, we employ 


Eq. j37l) 


and rewrite 



h(t) as an integral, in analogy to Eq. (f3"9"]l . This yields 
directly the magnetization (S z (t)) = h(t)/2, 



2D 



k{x) 



x 2 (l+jln(f)) 



{cos{xt) — \)dx , 
(42) 

with k(x) defined in Eq. 1(27)) . 

For not too low energies, the logarithmic correction 
j ln(x/D) can be neglected to leading order in j, and the 
coupling j can be considered as unrenormalized. In this 
way, the correct perturbative short-time limit of Eq. (I26[) 
is recovered. In the long-time limit, J °° dx cos(xt) . . . 

can be replaced by Jq^^ dx . . . , where f(t) is a func- 
tion with values of 0(1). Thus, the long-time tail of the 
magnetization curve is readily obtained as 



(S*(t)) = -[l/(lnt-l/j) + l+j + 0(j 2 )] 



(43) 
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FIG. 5: We show (S z (t)) for different isotropic couplings j. In 
panel a) the full numerical solution of the flow equations (full 
lines in both panels) is compared to a fit to the asymptotic 
analytical behaviour of Eq. (|43[) (dashed lines, fit parame- 
ters given in Table U]). Beyond the perturbative short-time 
regime, described by unrenormalized spin-flip scattering ac- 
cording to Eq. (|28[) . the logarithmic renormalization of the 
spin-flip scattering coupling J x sets in. This leads to a log- 
arithmically slow relaxation of the magnetization, and the 
asymptotic value {S z (t — * oo)) = 0.5 + | (dashed lines in 
panel b)) is reached extremely slowly. 



with the asymptotic value 

(S z (i^oc)} =0.5(1 + j + G(j 2 )) . (44) 

The steady state magnetization mirrors the behaviour 
observed in our toy model: the reduction from full polar- 
ization is j /2, which is twice the equilibrium value (see 
Eq.S5J. 

We conclude that Eq. (j42|) describes both the short- 
time and the long-time limit of (S z (t)) exactly. This 
result can be also derived within a more sophisticated 
ansatz for the flowing couplings J~'(i?), that parame- 
terizes also the flow of couplings considerably above the 
Fermi energy, as detailed in appendix lEl 



B. Anisotropic regime 

1. Flow of the exchange couplings 

Before turning to the dynamical behaviour of the ob- 
servable (S z (t)), we recall important scaling properties of 



anisotropic exchange couplings J x and j". In the ferro- 
magnetic regime, we have j" < and a stable fixed point 
of the scaling equations (fT0|) exists only if \J»\ > \ J ± \. 
In this case, we have 



j x = 0, (45) 

and according to Eq. (fT0|) . the dimensionless trans- 
verse coupling decays asymptotically as j ± (B) = aB^ / 2 , 
where a is a non-universal number. Numerical tests show 
that a is within good accuracy identical to the coupling 
j 1 - as long as |j" / > 2. These properties will be used 
in the following. 



2. Flow of the spin operator 

Our discussion of the flowing spin operator starts again 
with an analysis of the couplings jeg(B) that are gen- 
erated during the flow of the observable S z . As justified 
in appendix [Cj also for anisotropic couplings it is pos- 
sible to approximate the flow of the couplings jj;g(B) 
by Eq. (j37|) . In the limit of small energy differences 

Ae = |sg — eg | <C D, the integrated flow J °° Jf- R (B)dB 
is dominated by the slow asymptotic decay of J^ R (B), 
and we can set Jj^(B) = (a/ ' pf)B^ I 2 , yielding the ap- 
proximated couplings 



( £ k- £ k>) e 



o 

aAe 



N 



dB 



e ~BAs* B ^ dB 
Np F Jo 

asgn(Ae) 



N PF \Ae\^ 



+./" 



(46) 



Since Eq. (|38[) is valid also for anisotropic exchange cou- 
plings, we can again use it to calculate the renormal- 
ized coupling h, and by employing the couplings (|46[) in 
Eq. (|38|) . we obtain the result 



2L> 



dxa 



k{x) 

r 2+2jW 



1+ +G(J 2 ). (47) 
4? II 



3. Time-dependent magnetization 

We proceed as in the isotropic regime and plug the 
approximate couplings 7^^, from (|46p into the formal ex- 
pression from Eq. (|4"Tj) to obtain the coefficient h(t). This 
result is at least correct in the long-time limit. In this 
way, we obtain the time-dependent magnetization as 



(S'{t)) 



1 1 



2D 



a 2 ^L(l - cos(xt))dx. (48) 



r 2+2jll 



As argued before, for anisotropics </" /J x | > 2, it is pos- 
sible to set a fa j . In this case, again the short-time 
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FIG. 6: Short-time behaviour of the magnetization for 
anisotropic couplings from a fully numerical solution of the 
flow equations in comparison with the perturbative result 
from Eq. ([28} (dashed line). The perpendicular coupling is 
fixed to J x = —0.04, the curve with j" = —0.04 corresponds 
to the isotropic case. For increasing \J"\, the spin-flip scat- 
tering coupling J x is stronger renormalized, and the magne- 
tization tends to decay slower than predicted by the unrenor- 
malized perturbative result. 



TABLE II: For large times, we fitted the full numerical solu- 

tions of the flow equations against the function at V II +b, 
as depicted in Fig. [7K • For couplings — j" < 0.1, quantita- 
tive agreement with the analytical result (S z (t)) = 0.5(1 — 



7(2jU)^" + a 2 /(2]")+0(j 2 )) 



is very 



(S z (t)} = at V J il J i +b 


j 11 = -0.05 


j" = -0.1 


a 
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FIG. 7: Long-time result of the magnetization relaxation 
for anisotropic couplings and comparison with a long-time 
fit. All parameters are equivalent to those from Fig. [6] In 
panel a), the full numerical solution of the flow equations 
(full lines in both panels) is fitted against the analytical pow- 
erlaw behaviour stated in Eq. (|49p (dashed lines). At time 
scales where renormalization of J x by j" becomes signifi- 
cant and leads to a deviation from the perturbative result, 
again the asymptotic fit describes the relaxation process very 
well. In panel b), a comparison with the saturation value 
{S z {t -> oo)) = 1/2 + a 2 pF/(4j") (dashed lines) is given. 
Clearly visible, anisotropic couplings J x < j" lead to much 
faster saturation of the magnetization than isotropic cou- 
plings. Fitting parameters are given in Table ITT1 



limit of the flow equation result for (S z (t)) coincides 
with the exact perturbative result of Eq. (f25|) . Since 
our approximations become only exact at low energies, 
we have in this case no argument why they should in 
general hold in the short-time limit, where the influence 
of higher energy scales might invalidate our approxima- 
tions. In Fig. [6J the agreement of the numerical solution 
of the anisotropic flow equations with the perturbative 
short-time formula Eq. ([28]) is shown. Anisotropy due to 
increasing \J" \ does not influence the short-time limit of 
the observable (S z (t)), since j" enters this quantity only 
in third order of perturbation theory. 

The long-time tail of the magnetization curve can 
be obtained as in the isotropic case by replacing 
dx cos(xt) . . . with Jq^^ dx . . . , where fit) is a 
function with values of 0(1). This yields the asymptotic 



result 

<*■"»= K 1 -^'^ +o</) )- (49) 

In Fig. [3 the qualitative difference between the 
anisotropic and isotropic couplings is depicted. Clearly 
visible, anisotropy leads to qualitatively faster satura- 
tion of the magnetization to its steady state value. The 
algebraic long-time tail of the magnetization curve sets 
in immediately after the perturbative short-time limit. 
A numerical fit of the long-time tails yields very good 
agreement with the analytical prediction of Eq. (j49|) , as 
demonstrated in table [TTJ Furthermore, our calculations 
show that the steady state magnetization (S z (t — » oo)) = 

2 

is again reduced twice as much from full polarization 
than in equilibrium. 
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C. Conclusions 

The numerical results show that spin relaxation of 
ferromagnetic impurities is described by two different 
regimes, which are well described within the flow equa- 
tion method. An initial short-time regime, driven by 
unrenormalized spin-flip scattering with an amplitude 
oc J 1 - 2 can be described within unrenormalized perturba- 
tion theory. At low-energy scales (corresponding to large 
time scales), renormalization of spin-flip scattering slows 
down the relaxation rate. Quite naturally, isotropic cou- 
plings are renormalized only logarithmically in compari- 
son to algebraic renormalization of J 1 - in the anisotropic 
case, leading to much faster relaxation of anisotropically 
coupled spins. It is important to note that our results 
are to leading order in J independent on the spin quan- 
tum number S of the impurity spin. As detailed in ap- 
pendix [B] a generalization to arbitrary S leads to 

(S*(t)} = Sh(t), (50) 

where h(t) is the universal function which we calculated 
above. 



This factor two occurs in a general class of discrete 
systems^ in weak-coupling perturbation theory. In a 
simple exactly solvable toy model with two lead levels, we 
could explicitly evaluate this factor r(g) beyond its weak- 
coupling limit r(g = 0) = 2, see Fig. [2j The toy model 
results for positive (that is antiferromagnetic) coupling g 
are of course not relevant for the actual Kondo model due 
to the strong-coupling divergence that makes perturba- 
tion theory invalid for g > 0. On the other hand, the 
fact that the running coupling flows to zero on the ferro- 
magnetic side is just the reason why our analytical flow- 
equation results become asymptotically exact for small 
negative g. This is also supported by exact numerical re- 
sults using time-dependent NRG 4 which showed excellent 
agreement with the analytical flow-equation calculation 
for weak ferromagnetic coupling. One can even under- 
stand the differences between the exact numerical results 
and the perturbative flow-equation results in Ref. 0: Ac- 
cording to Fig. [2] one would expect the flow-equation cal- 
culation to overestimate the non-equilibrium reduction of 
the magnetization (since r(g) < 2 for g < in Fig. [2]), 
which is exactly what one observes in Ref. |3 for larger 
ferromagnetic coupling. 



VI. SUMMARY AND DISCUSSION 

We have discussed the real time evolution of the ferro- 
magnetic Kondo model initially prepared in the product 
state 



IV) = I t>® \FS), 



(51) 



where \FS) is the non- interacting Fermi sea. Using the 
flow-equation method, we have shown that the time- 
dependent magnetization (S z (t))^ approaches a nonvan- 
ishing asymptotic value (S z (t — ► oo))^ logarithmically 
(with a power law) in time if the model is isotropic 
(anisotropic). One key observation is the fact that 
(S z (t — > oo)}^ differs from the equilibrium value (S z ) eq 
for infinitesimal positive magnetic field. In particular, 
this implies that the system retains a memory of its ini- 
tial preparation for all times. 

While this might appear surprising at first sight, it just 
reflects the following well-known property of the equilib- 
rium model: If one prepares the system in the initial state 
([ST]) and then slowly switches on the ferromagnetic cou- 
pling to the leads, according to the adiabatic theorem the 
system will evolve to the equilibrium ground state with 
magnetization {S z ) Qq . However, if the initial state were 
| J.) ® \FS), then the asymptotic magnetization would 
be inverted. In this sense even the equilibrium model re- 
tains a memory of its initial preparation due to ergodicity 
breaking. 

In the weak-coupling limit we could prove that asymp- 
totic non-equilibrium magnetization vs. equilibrium mag- 
netization differ by a factor two: 



As discussed in Ref. 0, our results for the ferromag- 
netic Kondo model have direct applications to molecular 
quantum dots. For instance, the magnetic dynamics of 
single-molecule magnets coupled to metallic leads may 
be captured (in the co-tunneling regime) by a ferromag- 
netic Kondo models Another experimentally relevant 
situation are coupled semiconductor quantum dots with 
partially screened spins. 

On the methodological side, we have shown that the 
flow-equation method can be applied successfully to real- 
time evolution problems in quantum impurity models. 
We envision further applications of our method in the 
study of non-equilibrium dynamics near impurity quan- 
tum phase transitions^ 



(52) 
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APPENDIX A: MATRIX REPRESENTATION OF 
THE TOY MODEL 

The five basis states with total spin S = 1/2, S z = 1/2 
and two fermions in the electronic levels are given by 



l^o) = 


0} <g 


U>«? 


5 It) 


Ml) = 


T>® 


1 T>® 


II) 


\>h) = 


1) ® 


1 T>® 


It) 


IV>3> = 


U)c 


5 |0) € 


5 It) 


|^4> = 


T>® 


1 1) ® 


It) 



(Al) 



where the first ket in the tensor product corresponds to 
the c-electron in our toy model, the second ket corre- 
sponds to the ci-electron and the third ket corresponds to 
the impurity spin. The subspace spanned by these five 
states is denoted by H 2 i ■ Due to the SU (2)-symmetric 
interaction, the time-evolved initial state as well as the 
interacting ground state in the presence of an infinitesi- 
mal positive magnetic field lie in TL 2 i ■ 

In the subspace 7i 2 i , the Hamiltonian ([2]) of the toy 
model is represented by the matrix H = H$ + g Hi nt , 
H i+ltj+1 = (ifiilHtyj) for i,j = 0, 1,2,3,4: 



/ -2 \ 
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V 0/ 
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(A2) 
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(A3) 



The states \ips) and decouple, so that one effectively 
only has to diagonalize a 3 x 3-matrix to i) solve the dy- 
namics of the toy model exactly and to ii) determine the 
ground state magnetization. Since the diagonalization is 
a trivial step with lengthy expressions, we will not give 
details here. 

However, it is interesting to see explicitly how the fac- 
tor 2 in the ratio r(g) from |(7J) comes about in the weak- 
coupling limit. In perturbation theory the interacting 
eigenstates are given by 



(ipi\Hi nt \ipj) 
Ej — E.i 



+ 0(g 2 ), (A4) 



where Eg = —2, E\ = 2, E% = are just the eigenval- 
ues of Hq. For the observable O = S z — 1/2 one has 
(-00 1 l^o) = (Tpi\0\ifti) = and therefore the ground 



state expectation value 

o eq = (<>Po\d\i> ) 

= 5 2 (^2|0|^2) 



\(^2\H mt \^ )\ 2 
(E 2 - So) 2 



0(g 3 ) 



The time-averaged expectation value starting from the 
initial state \ipo) is according to © 



(A5) 



and one can easily verify that only the terms j = 0, 2 con- 
tribute plus corrections in order g 3 . Each of these terms is 
identical to O eq plus again corrections in order g 3 , which 
proves 



r(g) = 



0(t) 



a 



= 2 + 0(g) 



(A6) 



i-q 



The general proof in Ref. la is a generalization of this 
argument. The universal factor 2 in the weak-coupling 
limit plays e.g. a key role in the thermalization of a Fermi 
liquid after an interaction quench, that is the opposite 
limit of the adiabatic Landau Fermi-liquid paradigmi^ 



APPENDIX B: FLOW EQUATIONS FOR 
GENERAL SPIN S 

We now analyze how the magnetization curve (S z (t)} 
is modified if the impurity spin has a general quantum 
number S > h- The commutation relation [S i , S-*]- — 
itijkS k is independent of S, and we will show that also 
the flow equations for the spin components S a remain 
unchanged to leading order in J. We illustrate the proof 
only for the case of isotropic couplings, since the prove for 
anisotropic couplings requires only slight modifications. 

Using the ansatz (j2"D)) together with the generator (fl"5|) , 
it is readily seen that modifications to the flow equations 
for S a (B) can only arise due to the commutator 



[:S-Sfp :,(:Sx s,-,. :) a 

l tv ' >• u'u > 

We now decompose 

S $ — ^ij 



k 



(Bl) 



(B2) 



the operator obtained by project- 
^ 6ijk.S k as a basis operator 



where %Y,k e vkS k is --■ 

ing .S"'.S'' on the operator Y^ k , 

in spin Hilbert space and Xij is the operator defined by 
Eq. (|B2p . Using this decomposition in Eq. (|B1[) . we ob- 
tain 



[^•%:,(:Sx^:n_ = 



1 i 

+ 42 



Y, a U (X >-^'S k '[. ct a cl p :,: ct c tt „ :]+) . 
k' ' 

(B3) 
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It is now readily seen that corrections to the flow equa- 
tions Eq. (|2"Tj) can only arise from the term in Eq. (|B3[) 
that contains the operator Xij . This term can be rewrit- 
ten as 



ijkl 



(B4) 



and due to the decomposition assumption (|B2[) , it is read- 
ily seen that this expression has a vanishing projection 
on the operator (: S x s~, - :) a . Thus it can only be gen- 
erated as a subleading correction in 0(J 2 ) of the flowing 
spin operators. From the identity J^. xn — S(S + 1), it 
is finally seen that the flowing spin components S a have 
corrections which are of 0{ J 2 S 2 ). For the magnetization 
curve (S z (t)), this means that the relation 



(S'(t)) = 2S(S'(t))\s=l+0(J*S J ) 



(B5) 



is fulfilled. 



APPENDIX C: VALIDITY OF THE 
TREE-LEVEL APPROXIMATION 

In section V, we made use of the approximative differ- 
ential equation 

^§1 =h{B)[e 9 -e n )J^{B)+0{{e 9 -e % )J^) . (CI) 

In case of anisotropic couplings \ J ± \ < \ J"\ and j" < 0, 
the couplings J- - contained in the neglected terms flow 

k'k 

to a finite value at the Fermi surface (given by e^ = 
e£,=0), while Jg^B) will flow to zero. A justification 

of neglecting all terms of 0(J 2 ) in 
consider the flow equations including all terms of up to 

o(J 2 ), 



1 1 ; " dlk ' k has therefore to 



dB 



N 



1 



D 



x e~ 



de sgn(e) 



(Ae + eo - e 



N 



S(Ae+££ '" e)2 / dB'(e-e % ,)Ji- R (B') 
Jo 



{ep - e) ^l e -B(e-e p r 



N 



x 



[ dB'(Ae + e p -e) 
Jo 

jj- R {B')e- B '^^-^ 2 



(C2) 



where we defined Ae = e? — er, and made use of the in- 

k k 



frared parametrization (|32p . Only in three integration in- 
tervalls of width Ae, the integrand in Eq. (|C2[) is nonzero 



and we can simplify according to 



d 1k'k = A c J m( B ) -B(Ae) 2 

N 

[\B\e-e % ,)ji R {B')e- B '^ 2 ] 
Jo 



dB N 
1 

+ 2 



+ 0(er BD *) , 



(Ae+ ££ ,- e ) 2 



(C3) 



where 0(e BD ) arises from integrals close to the band 
edges. Since the perpendicular coupling Jj~ R (B) asymp- 
totically decays as Jj~ R (B) oc j^B^/ 2 (see sec. IV B|) . a 
lower boundary for the contribution of 0( J) in (|C3j) is 



PF 



(C4) 



while it is readily seen that the contribution of 0(J 2 ) has 
the upper boundary 



BJ x jWO(Aef 



(C5) 



as long as B < (Ae) 2 . Hence, in the regime B < (Ae)" 



we can approximate d ^ k R k ~ AeJ- L (B)e B ( Ae ) if 



Ae 2 B 1 ~ j /2 < - 
j 



1 



(C6) 



The remaining regime B > (Ae) 2 is described by ex- 
ponential decay of 3b* an( i i s n °t of importance. In 
conclusion, we showed that the approximation 

is valid also for anisotropic couplings with \ J ± \ < |j"| 
and JH < 0. 



APPENDIX D: AN ALTERNATIVE ROUTE OF 
CALCULATION FOR EQ. (|23)) 

An alternative derivation of the formal result (|4"Tj) for 
the magnetization (S z (t)) — h(t)/2 exploits the sum 
rule ([22]) . Since the operators S z (t) and S z (t) are re- 
lated to the flowing operator S Z (B) by a unitary trans- 
formation, the sum rule (S Z (B) 2 ) = 1/4 applies for both 
operators, and by subtracting the corresponding sum rule 
for expression (|30|) from that of expression (|3ip . the co- 
efficient hit) can be expressed as 

h 2 it) = h 2 +j2(\%~,m 2 -h*i(t)\ 2 ) 



kk' 



x n(fc')(l - n(fc)) + C(J 2 ) . 



(Dl) 
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Eq. (jDip is another starting point to calculate the mag- 
netization h(t)/2. The coefficients 7g/g(u and 7g/g(i) 
in (|Dlj) are given by Eqs (|30|) and (|40| , and we obtain 



fr 2 (t) = ^ 2 + 5>gg-,| 2 (-l + 2cos[( £ g-eg-,)t 
x n(k')(l -n(k)) + 0(J 2 ) . 



(D2) 



Using the known result h — 1 + O(J), it is easily seen 
that the right hand side of Eq. (|D2[) behaves as 1 + (J), 
and the expansion of ^Jh 2 {t) to first order in the coupling 
J is therefore 

fc(t) = ^ + ^l7^| 2 (cos[( ££ -£, r ,)t] - - 



x n(fc')(l -n(fc)) + C(J 2 ) 
This result is equivalent to Eq. Ell) . 



(D3) 



APPENDIX E: DIAGONAL 
PARAMETERIZATION OF ISOTROPIC 
COUPLINGS 

In section |V1 we integrated the flow equation 



dB 



(El) 



by using the approximative parameterization ji., (B) 



(J IR (B)/N)e- 



This parameterization assumes 



isotropic couplings J^{B) = jiL (£?) and requires low 
energies |eg| <C D , \ep | <C D. Here, we provide a more 
general integration of Eq. (IEip by parameterizing the flow 
of the couplings Jtg (B) by an ansatz that also describes 
couplings above the low-energy limit appropriately. One 
makes the ansatz 



Jkk'(B)=J W (B 



-B(e s - 



(E2) 



with Jjjj(B = 0) = J/N < 0, where k'k is a label for the 
energy median 



k'k 



(E3) 



such that the couplings J-g~{B) depend only on the en- 
ergy e-gp- and B. This approach has been dubbed diag- 
onal parameterization, since it approximates the full set 
of N x N couplings by its N diagonal entries. It has been 
used previously in the context of the non-equilibrium 
Kondo model in Ref. [l|| and a detailed discussion is 
given m Ref. [H It is trivial that this ansatz becomes 
asymptotically exact for the low-energy couplings near 
the Fermi surface, since the coupling Jj=y(B) will reduce 
to the IR-parameterization Jm(B) used in sec. Mil once 
e k ' e k> 0- In addition, it describes the deviations from 



the flow of Jm(B) caused by interactions between matrix 
elements of higher energy scales. 

The number of flow equations for the couplings reduced 
now from N x N to the N differential equations (using 



J k ~ J kk> 



dJ^B) 
dB 



^2 (e , + £ ,)J^e- 2 ^-^K / T)-i] 



0(J 2 ) 



(E4) 



Off-diagonal couplings with huge energy differences de- 
cay exponentially fast, and it is possible to make the 



approximation J-^r(B) 



Jg(B) in Eq. (|E4)l . Further- 
more, the initial condition for Jg(-B) can be adjusted to 
J^(B = D~ 2 ) = jr, since the flow of Jj:{B) is unrenor- 
malized for energies near the band cut-off, corresponding 
to B < D~ 2 . Now, the diagonal couplings flow approxi- 
mately as 

dJ%{B) 
dB 



;) 2 JW--(n(/i) - -) 

kp pk y VPV 2 



J f^V 2 -t 



(E5) 



2BN 

where we evaluated the momentum sum at T = 0. This 
differential equation can be integrated in the closed form 

-l 



J 



PF 



B 



dB 



/ exp(-2B'e| 



D- 



B' 



(E6) 

In order to evaluate the couplings 7gg; by integrat- 
ing Eq. (|E1|) . we approximate the couplings </g(-B) as 
Jg(-B = (eg — eg)~ 2 ). This is justified up to corrections 
of 0(J 2 ), since J^{B) depends logarithmically or weaker 
on the parameter B. We can further evaluate the cou- 
plings J^(B = (eg — £g)~ 2 ) by performing the sequence 
of approximations 



dB 



exp(-2B'el) 



B' 



dB 



,exp(-2i3'( eE - e p) 2 ) 



D- 2 
oo 



B' 



Oil) 



dB 



,exp(-2B'(ej 



d- 2 



In 



B 



D 2 



) + o( 



2fe 



0(1) 



D 2 



0(1) . 
(E7) 



We assumed in the second line that eg and eg — eg, are of 
the same order, since the phase space where this is not the 
case has subleading measure in momentum sums gg ; . 



Assuming furthermore In 



D 2 



> 1 we made use 



of the asymptotic expansion of the exponential integral 

/•oo t 

Ei(l,x)= —dt = - 1 -\nx + 0{x) 
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in the last line, with 7 « 0.5772. 

The parameterization (|E2[) finally justifies that the co- 
efficients 7£ fe -, can be obtained using the infrared parame- 
terization J4, (B) = {J IR {B)/N)e- B{ - € k-H') 2 in the limit 

|ln ^ £fc ) 1^1- This result is of particular impor- 
tance for the the time dependence of the magnetization 
(S z (t)), which is fully described by the couplings 7^ 
(see Eq. gl])). 



APPENDIX F: NORMAL ORDERING 

Within the flow equation approach, fermionic oper- 
ators occuring in transformed observables are usually 
normal-ordered in order to separate the various inter- 
action terms generated during the flow into irreducible 
objects in the spirit of a diagrammatic expansion using 
Wick's theorem.— In practice this procedure is imple- 
mented with respect to a given reference state or density 
matrix, such that normal ordering of a product ct c r , 

k(T k' (7 

of fermionic operators amounts to subtracting the expec- 
tation value with respect to the given reference state or 
density matrix, 



c ka c k'. 



ct c-, , 

ka k a 



(ct (U ) 

x ka k'a'' 



(Fl) 



where we introduced a general multiindex I comprising 
both spin and momentum label. Since a calculation of 
the full dependence of nw (B, t) on the parameters B and 
t turns out to be a very complicated problem of its own, 
we show here that the improvement achieved by includ- 
ing this dependence in the contractions nw (B = 0, t = 0) 
enters the flow equations for spin operators (the Hamil- 
tonian) only to 0(J 2 ) (to C(J 3 ) ), having no influence 
on the leading order of our calculation, which is 0(J) 

(o(J 2 )). 

In order to reduce (|F3|) to the form nw{B,t) = 
nw(B = 0,t = 0) + O(J), we write the transformed 
operators c t (B) = W(B)ciU(B) in the general form 

c, (B) = //(S)c, + J x composite operator+£>( J 2 ) , (F4) 
with a flowing coupling parameter fi(B) and the unitary 
transfomation U(B) generated by the generator (|15p . 
Here and in the following we exploit that fi(B) behaves 
as 1 + O(J), since ff{B) determines the quasiparticle 
weight of the conduction electrons that is renormalizcd 
only perturbatively in J. In this way the contractions 
nu'(B,t) simplify to 

n w (B,t) = (^\e tHt cl Cll e- iHt \^) + 0(J) . (F5) 



In our particular problem, we used the initial state \ipi) as 
reference state for the normal ordering procedure. How- 
ever, this reference state will look very different at some 
finite time t or if it is transformed into the representa- 
tion of the flowing operator by the unitary transforma- 
tion U(B) induced by the flow equation approach. In 
order to use a normal-ordering prescription with respect 
to the same quantum state for all values of B, we need to 
define normal ordering with respect to the time-evolved 
initial state in the basis for a given value of B, 



iHl 



del 



mB,t)) = U(B)e mt \^) 



(F2) 



The contraction of fermionic operators then depends 
both on B and on time t, 



n w (B,t) d = {ij>i\e 



-iHt 



W(B)c] Cl ,U(B)e lHt \^) , (F3) 



can 



The time-evolved fermionic operator e c, e 
now be formally calculated in the forward-backwards 
transformation scheme depicted in Fig., generating an ef- 
fective solution of the Heiscnberg equation of motion for 
c ; of the form c ; (t) = di(£)c ; + J x composite operator + 
0(J 2 ), with a coupling function di(t) = 1 + O(J). Now 
the contractions can be written as 



nu>(B,t)(ip(B,t)\clc ll(r ,\ip(B,t)) = (^|cJc;,|^) + 0(J) . 

(F6) 

Since nw(B,t) enters our flow equations only in terms 
of highest considered order in J the correction of 0( J) 
occuring in (|F6[) enters our calculation only to subleading 
power in J. Therefore, it is sufficient in our calculation 
to use the simplified contractions nw — {ipilcj c v for 
normal ordering. 
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